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Abstract 

Modulational dynamics of oscillatory travelling rolls in magnetoconvec- 
tion is studied near the onset of a Hopf bifurcation. Using weakly nonlin- 
ear analysis, we derive an envelope equation of oscillatory travelling rolls in 
. the plane perpendicular to an ambient vertical magnetic field. The envelope 

equation is the Davey-Stewartson (DS) equation with complex coefficients, 
from which we obtain criteria for the modulational (Benjamin- Feir) insta- 
q ■ bility of oscillatory travelling rolls. 

czj ' 

1 Introduction 

Or 

A variety of spatially and temporally periodic patterns is found in weakly non- 
^ ■ linear Boussinesq convection in an imposed vertical magnetic field [|I|]. However, 

nonlinear modulation of the periodic patterns in the horizontal plane has thus far 
received less attention. Two kinds of bifurcations are known to convective pat- 
terns: one is to steady patterns and the other is to oscillatory patterns (a Hopf 
bifurcation) [[![]. Near the onset of the former bifurcation, modulational dynam- 
ics of travelling rolls may be described by the Newell- Whitehead-Segel (NWS) 
equation in terms of the Newell- Whitehead (NW) orderings ||. Near the onset 
of a Hopf bifurcation, the same ordings may yield a more complicated envelope 
equation H, which includes various terms of different orders. 
In this paper, introducing consistent orderings (different from the NW orderings) 
to an envelope equation near the onset of a Hopf bifurcation, we derive the DS 
equation with complex coefficients, where horizontal incompressible flows couple 
to oscillatory travelling rolls in magnetoconvection. The analysis of a instability of 
its spatially uniform solution yields criteria for the modulational (Benjamin- Feir) 
instability of oscillatory travelling rolls. 
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2 Hopf Bifurcation 



Boussinesq convection in an imposed vertical magnetic field is described by the 
equations QTJ| 

d t u + (u-V)u = -VP + aRQz + a(Q(V x B) x B + aAu, (2.1) 
d t Q + (u-V)& = w + AQ, (2.2) 
d t B + (u-V)B = (B -V)u + (AB, (2.3) 

together with 

V ■ u = V • B = 0. (2.4) 

Here n = (u, v, w) is the dimensionless velocity in (x,y,z) coordinates, z is the 
unit vector in the vertical (z) direction, is the dimensionless temperature de- 
viation from the conduction state and B is the dimensionless magnetic field. 
The parameters are the Prandtl number a, the Rayleigh number R, the Chan- 
drasekhar number Q and the ratio of ohmic to thermal diffusivity denoted by £. 
The magnetic field is given by 

B = z + b, 

where b = (a, b, c). The boundary conditions are 

d z u = d z v = w = Q = a = b = at z — 0, 1. 

The linear stability analysis of the conduction state u = G = b = shows that 
oscillatory convection sets in, for ( < 1, at 

P o (vr 2 + P) 3 + 0(1 + , vr 2 + fc 2 C(^ + C) 2p 

Here k is the horizontal wavenumber determined by minimizing the critical Rayleigh 
number R . Thus k = k , where 

(7r 2 + ,2X3 _ 3 2( 2 , 2)2 = ^ 4 

V 0J 2 V o; 2(l + cr)(l + 



3 Envelope Equation 

In this section, we derive the equation which describes the nonlinear evolution of 
a slowly varying envelope of oscillatory travelling rolls near the critical Rayleigh 
number R . For values of R close to i?o, 

R = Ro + e 2 R {2) (e<l), 



2 



we investigate the weakly nonlinear evolution of the wavepacket centered at the 
critical wavenumber k and the corresponding frequency u = u(k ) by approxi- 
mating u, b, and P as 



eb 1 ^,r ] ,r)E 



/ COs(7Tz) 

cos(7T2;) 
y sin(7T2;) 

/ sin(7rz) 



(c.c.) + eV 2 ) + eV 3 ) + 



(c.c.) + e 2 & (2) + e 3 b (3) + 



(3.1) 



(3.2) 



sm{7(z / 

\ C0S(7TZ 

6 = eQ 1 {£,r ] ,T)Esm{7rz) + (c.c.) + e 2 9 (2) + e 3 6 (3) + •■■, 
p = P 1 (^ V} T)Ecos(7Tz) + (c.c.)+e 2 P {2) + e 3 P (3) + •••, 
E = exp[i(koX — Uot)], 

where (c.c.) denotes the complex conjugate of the previous term, and 

£ = e (x - At), rj = ey, r = e 2 t. 

The present ordering ( |3.5| ) of the stretched variables £, rj, r is the same as one 
introduced in the derivation of the DS equation ||. Substituting the expansions 
nj)-(|3.4j) into EqsfETlT)- (12. 4|) and using the method of multiple scales H, we 



(3.3) 
(3.4) 



(3.5) 



into Eqs( |2.1| )- ( |2.4| ) and using the method of multiple scales 
obtain from the leading order equations(linearized equations) 

-2 ■ , \-l 



ci 
Pi 

Ui 



L e w 1 , L e -■ 
irL c wx, L c 
irLpWi, L t 



y K — IUq 
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i-l 



-a(R L e /K 2 + CQL C ) 
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l—Wi, 

k 

vi = 0, 
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a t = -i—c\, 
k 



where k 2 = k 2 + n 2 and the linear dispersion relation becomes 



f = <7K — iu — Ti L p — aR L e = 0. 
The second order field variables are expressed as follows 



9 (2) 

c< 2 > 

p(2) 



u 



(2) 
,(2) 



Esm.(irz), 
e {2) Esm(7rz) + (c.c.) + sin(27rz), 
[cS 2) cos(Trz) + c { 2 ] E]E + (c.c), 
[p[ 2) cos(tt2) + {P 2 (2) cos(27T^) + P 2 (2) }P]P ^ 

+P (2) COs(2tTz) +P (2) , 

Ui E cos(ttz) + (c.c.) + Uq cos(27rz) + v!q\ 
di E sm(7rz) + (c.c.) + sin(27T2;), 



Vi E cos(jrz) 



v (2 \ = b?Esm(nz) 



c.c, 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 
(3.11) 
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where and v$ are determined by 



7T 
7T 2 



'l-aCQ^L^d^Wr 



(3.12) 



The other second order amplitudes such as 0i , \ and so on are given in 
Appendix and give 

A = [d k u(k)] k=ko = d k uj . 

As shown in Appendix, the solvability condition for the third order variable w\ 
yields the following equation of the envelope of the first order vertical fluid velocity 
wi : 

2 * 2k v \d w J J d w J 

(3-13) 

where d k u is real, while S^Uq = (dlu) k=ko has a complex value in general and 



q = -Ro^ 2 LeL' e + aQ^—- 

SiaQir 6 
k\Q + A^) c c 



3k 2 - 7T 2 



\L C \ 2 - Lc 2 



iQn 2 



Q + 47T 2 



7i 2 a(QL 2 c + aRo-^Ll + 



k 2 - 3tt 2 



K 



(3.14) 



A coupled system of equations ( |3.13|) and ( |3.12| ) is the DS equation with complex 

_(2) 

coefficients, in which the real field u represents z independent incompressive 
horizontal flows varying slowly. 



4 Modulational Instability of Oscillatory Trav- 
elling Rolls 

For later conveniences, the DS equation with complex coefficients (a coupled 
system of equations ( 3.13| ) and ( |3.12 ) ) is rewritten as follows. 



i dt%> + ad 2 x m + pd 2 ^ + (7|^| 2 + su)^ = ir^, 
(d 2 x + ad 2 )u = (bd 2 x + cd 2 m 2 . 



(4.1) 
(4.2) 
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where a,/3 and 7 are complex constants (a = a' + ia",/3 = + 7 = j' + ij"), 
while s, r, a, b and c are real constants (a = 1, b = 0, (5" — in the present case). 
A spatially uniform solution of Eqs.(|4.1|) and (|4.2|) is given by 



* = ^0 = V'o exp(-ifit), w = 0, 

where |?/>o| 2 = r/j" and = — 7'|^ | 2 - Setting 

* = ^0 + Vi(*)exp[i(fe • r - n x t)] + ip 2 (t) exp{i(-k ■ r - tt 2 t)), 
u = Ui(t) exp[i{k ■ r — (Q 1 — fl)t}] 

+u 2 (t) exp[i{-k ■ r - (0 2 - Q)t}] + (c.c), 

where A; = (k x ,k y ),r = (x,y),2Q = ^1 + ^2 and linearizing Eqs.( f4.1|) and (|4.2| ) 
with respect to ipi and ip 2 , we have 



= [r-*{c^ + ^-ni-(7 + ^)|^rM + <7^^, (4-3) 
[r + z{a*fc 2 + - n 2 - (7* + 7*)l^o| 2 M ~ »7>S Vi, (4-4) 



^2 _ r,. ;f ,, : /.2 ?i /.2 n .1 \ 1 . , 1 2 VI , , - i , ,.i.2 



where 7 = 7 + s(6A^ + cky)/(k 2 + ak 2 ) and * denotes the complex conjugate. 
The linear equations ( |4.3| ) and ( |4.4j ) have exponentially growing solutions if the 
following condition is satisfied. 

{a'i + a'Y)k 2 x + (/?Y + (3"i')k 2 y > 0. (4.5) 

Since a > in the present case, Eq. (|4.5|) yields the following instability criteria. 

a = 0/(7' + sb) + a"j" > 0, (4.6) 

or 

$ = + sc/a) + £"7" > 0, (4.7) 

or 

(q - (3/a) 2 + s[2(d + /3/a) + s] > 0, and d + /3/a + s > 0, (4.8) 

where s = s(c/a — b) (a' — /?'/a) > 0. The criterion ( (4.6| ) or ( J4.7|) is essentially 
the same as the case of the two-dimensional complex Ginzburg-Landau equation 
(Eq. ( |4.1| ) with s = 0). A new criterion ( |4.8j ) comes from the coupling between 
convetive rolls and horizontal imcompressible flows. If a < 0, although this is not 
the case in magnetoconvection, a spatially uniform solution of Eqs. ( |4.1| ) and 
(|4.2|) is shown to be always modulatinally unstable. 
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5 Concluding Remarks 

In this paper, we have derived a envelope equation of oscillatory travelling rolls 
near a Hopf bifurcation, which is not a type of the NWS equation || but the 
DS equation with complex coefficients. The NWS type equation has a defect if 
the group velocity at the critical (bifurcation) point does not vanish. That is, 
it consists of different order terms: the main term is linear and proportional to 
the group velocity, while the other interesting terms such as nonlinear terms are 
of higher order. Although our derivation is based on weakly nonlinear analysis 
with multiple scales similar to || , the present ordering of stretched variables (|3.5| ) 
is different from the NW ordering and yields the DS equation with complex 
coefficients which consists of the same order terms. 

Analyzing the modulational instability of a spatially independent oscillatory so- 
lution of the DS equation with complex coefficients, we obtain criteria of the 
modulational instability, which include not only the known criterion for the com- 
plex Ginzburg-Landau equation ,which was first given in ||, but also a new 
criterion due to the coupling between convective rolls and horizontal incompress- 
ible flows. 



Appendix 

A Higher Order Amplitudes 

(2) 

In terms of the first order amplitudes and w\ , the second order amplitudes 
defined in Eqs.( [3.6|) - (|3.11|) are given by 



r\( 2 ) T ( 2 ) T a (2) r (2) -fa 

= Lgw{ — tLeO^Wi, c\ = ttL c w\ — t7rL c o^Wi } 

'Ml 'M) 

Kq Kq Kq Kq 

Lq = \8^ Lq + dk Lg, etc., 



and 



(2) - L ' 9 \vn\ 2 7/ (2) - W * r"\v-\ 2 n {2) - 27l v {2) 

D (2) , aR V L 'e,\ |2 , °"C<5 k2 , ,2 

p ° = (i + -^r)M +^-N > 

5(2) _ Aj7 (a)_l!u. .2_fCQ ri -^U| 2 
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°CQk 2 2 5(2) ^ 7T K 2 7T 2 2 



^(2) _ 

2 ~~ o n ,2 
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C, 



2(k^ — iuo 



wiCi, Lc" = Im(L c ), L' g = Re(Le 



Lengthy but straightforward calculations of the third order terms in Eqs. (|2.lQ - 
( |2. 1| ) give the following third order field variables (envelopes) proportional to E 
in terms of W\ and 

ei 3) = L e w{ 3) -it e d^ ] -(Le^dlw.-L^-d^ 

-ik Ll(u { Q } - u£ ] /2)wi + TcLgQ^wx, 
c? = 7r[L c «;[ 3) -i4M 2) - (4/2)^«;i -^(a.-C^i] 

+ihL c [a£ ] /2 - 7tL c (m[, 2) + uf ] /2)}wi - 7rL c c[, 2) W-i, 
P{ 3) = 7t[L p w? - zL p 9^i 2) - (V2)cf Wl - l/(2fc )^ L p 5> 1 

+a(Qn[L c - L*(k 2 - 3n 2 )/ k 2 ]^ w_i 
-{(iaCQ/2){(k 4 + 4k 2 7i 2 - 7T 4 )/kK 2 }L c a^ + k^d^L^ 
-(tk n/2){a(QL 2 c - (aR /n 2 )L 2 e - (A/n 2 )}u^ 



+(aR n 2 /K 2 )Lo&Z ) ]w 1 , 



and 



+ (fco^ /4 2) - d R JR^)wt - q\ Wl \ 2 w u (A.l) 
where g is given in Eq. ( |3.14| ). Since 

/ = / = o, §4 = _L = -^v, 

Eq.pn|) gives Eq.flOf). 
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